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ABSTRACT 

We here consider the pressure equilibrium during an intense magnetic flare above the 
surface of a cold accretion disk. Under the assumption that the heating source for the 
plasma trapped within the flaring region is an influx of energy transported inwards with a 
group velocity close to c, e.g., by magnetohydrodynamic waves, this pressure equilibrium 
can constrain the Thomson optical depth tt to be of order unity. We suggest that this 
may be the reason why tt ~ 1 in Seyfert Galaxies. We also consider whether current 
data can distinguish between the spectrum produced by a single X-ray emitting region 
with tt ~ 1 and that formed by many different flares spanning a range of Ty. We 
find that the current observations do not yet have the required energy resolution to 
permit such a differentiation. Thus, it is possible that the entire X-ray/7-ray spectrum 
of Seyfert Galaxies is produced by many independent magnetic flares with an optical 
depth 0.5 < r T < 2. 

Subject headings: acceleration of particles — black hole physics — magnetic fields - 
plasmas — radiation mechanisms: non-thermal — galaxies: Seyfert 
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1. Introduction 

Although the X-ray emission from Seyfert galaxies 
is still not fully understood, there are indications that 
the mystery is starting to unfold. It is now believed 
that the emission mechanism is most likely thermal 
(see, e.g., Svensson 1996; Zdziarski et al. 1996). A 
flattening of the spectrum above about 10 keV (e.g., 
Nandra & Pounds 1994) strongly suggests the pres- 
ence of a component due to the reflection of the in- 
trinsic X-ray spectrum by nearby cold matter, which 
means that the X-ray emitting region is probably 
located above the accretion disk. In addition, the 
fact that the UV and X-ray luminosities from these 
sources are comparable and that their time variations 
are tightly correlated suggests that the UV and X- 
ray emitting regions are located very close to each 
other. There are also indications that the broad and 
redshifted Fe Ka line originates within the cold ac- 
cretion disk which extends inwards to the last stable 
orbit (e.g., Iwasawa et al. 1996). 

Many of these observational characteristics fit within 
the currently popular two-phase patchy accretion disk- 
corona model (Haardt & Maraschi 1991; Haardt & 
Maraschi 1993; Haardt et al. 1994), in which the X- 
ray emitting region is located within 'active regions' 
(AR) above the disk, where the energy dissipation 
rate greatly exceeds the local gravitational dissipa- 
tion inside the disk — at least for the lifetime of the 
region. Recent calculations by Zdziarski et al. (1996) 
show that the Thomson optical depth tt in Seyfert 
galaxies must be close to unity. From a theoretical 
perspective, this is highly significant given that the 
sample of sources covers a broad range of black hole 
masses, accretion rates, and other system parameters. 
Thus, it is not easy to understand why tt ~ 1, espe- 
cially since the energy transport mechanism has not 
yet been identified. 

Haardt et al. (1994) suggested that the ARs 
may be magnetic flares occurring above the accretion 
disk's atmosphere and showed that their compactness 
I may be quite high (~ 30), so that pairs can be cre- 
ated. The local compactness parameter I is defined as 
I = Lot I ' Rm e c 2 , where L is the source luminosity, R 
is the source size and ot is the Thomson cross section. 
In principle, it is possible to obtain still larger values 
for the compactness parameter (I ~ few hundred; see 
Nayakshin & Melia 1997a), thus creating enough pairs 
to account for the observed tt — l.Olg'2 (Zdziarski 
et al. 1996). This explanation for the observed value 



of tt ~ 1 based on the pair equilibrium condition, 
however, relies on the assumption that the particles 
are confined to a rigid box, so that no pressure con- 
straints need to be imposed. This is unphysical for 
a magnetic flare where the particles are free to move 
along the magnetic field lines. Therefore, as far as 
the two-phase model without a proper pressure equi- 
librium condition is concerned, the Thomson optical 
depth is a parameter, rather than a calculable quan- 
tity. 

In this Letter, we attempt to address this question 
of pressure equilibrium during an energetic magnetic 
flare above the surface of a cold accretion disk. We 
show that under the assumptions (i) that the energy 
is supplied to the X-ray emitting region by magne- 
tohydrodynamic waves with a group velocity close to 
the speed of light and (ii) that the emitting region has 
a very large compactness parameter I ^> 1, the pres- 
sure equilibrium between the outwardly directed X- 
ray flux and the inwardly directed MHD waves leads 
to the requirement that tt ~ 1-2, in accordance with 
the apparently "universal" value of tt seen in Seyfert 
galaxies. We will also show that both assumptions are 
likely to be valid during a physically plausible mag- 
netic flare above the surface of a cold accretion disk. 

2. The Energy Supply Mechanism 

The 'universal' X-ray spectral index of Seyfert 
Galaxies (e.g., Nandra & Pounds 1994) suggests that 
the emission mechanism is thermal Comptonization 
with a y-parameter close to one (Haardt & Maraschi 
1991; Haardt & Maraschi 1993; Fabian 1994). For 
the emission to be dominated by Comptonization, the 
compactness parameter needs to be large (I 3> 1; e.g., 
Fabian 1994). Observations of Seyfert Galaxies point 
to a global compactness parameter ~ 1 — 100 (Svens- 
son 1996 and references cited therein). 

For electrons and protons at a single temperature 
T e and an electron number density n e with the as- 
sumption of neutrality, the gas pressure is 2n e kT. 
The radiation energy density may be recast in terms 
of the luminosity L of the source, and therefore its 
compactness I. The total pressure is 
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P=^[2r T Q e + l(l + r T )/3}, (1) 

where 9 e = kT e /m e c 2 is the dimensionless electron 
temperature. The Compton y-parameter for thermal 
electrons is y = AQ £ tt{1 + 7t)(1 + 46 e ) (see Rybicki 



2 



& Lightman 1979) and is of order 1. Thus, since the 
dimcnsionlcss gas pressure in Equation (1) is always 
smaller than y, the radiation pressure dominates over 
the gas pressure in a one-temperature plasma when 
/> 1. 

One consequence of this is that the amount of en- 
ergy escaping from the source even during one light 
crossing time is larger than the total particle thermal 
energy. Thus, there must be an agent that energizes 
the particles to enable them to radiate at this high 
rate, and the presence of this agent must be dynam- 
ically consistent with the state of the system. We 
foresee two possibilities for the nature of this 'agent': 
(i) the gravitational field, and (ii) an external flow 
of energy into the system. These two cases are quite 
distinct physically. 

Insofar as the first possibility is concerned, the 
gravitational potential energy of the plasma (primar- 
ily that of the protons) is dissipated as the gas sinks 
deeper into the well of the black hole. The gravita- 
tional field does not provide a pressure, but it does 
compress the gas. However, this leads to an internal 
(radiation plus gas) pressure that varies from source 
to source as the physical conditions change. There 
does not appear to be a scale that sets tt to have a 
value of 1. For example, in standard accretion disk 
theory, the inner radiation pressure-dominated regime 
has an optical depth that depends on several param- 
eters, such as the accretion rate and the a-parameter 
(Shakura & Sunyaev 1973). The a-parameter reflects 
the rate at which the protons 'use up' their poten- 
tial energy, and so a change in this rate leads to a 
change in the equilibrium optical depth. It is even 
less obvious why tt should be <~ 1 in the gas pres- 
sure dominated regimes since there the pressure has 
no reference to tt at all. It seems that when the 
pressure equilibrium is dictated by the gravitational 
field (e.g., due to a compression of the X-ray emit- 
ting region), the Thomson optical depth should span 
a range of values depending on the source geometry, 
the specific parameter values and the particle inter- 
actions assumed to operate in the source. 

This is not so when the energy is supplied to the 
X-ray emitting region by an inflow of energy, e.g., via 
a magnetic field. The principal difference between 
the two cases is that the dynamic portion of the mag- 
netic field supplies a "ram" pressure that is related in 
a known way to its energy density. If the magnetic 
energy flux into the X-ray emitting region is known, 
this also constrains the inwardly directed momentum 



flux (the compressional force) into the system. Thus, 
the compressional force exerted on the active region 
by the magnetic field is expected to correlate with the 
source luminosity. What makes this useful in terms of 
setting the optical depth of the system is that a simi- 
lar correlation exists between the luminosity and the 
outwardly directed radiation pressure in the emitting 
region. But in this case, the pressure also depends on 
tt- Assuming a spherical geometry for simplicity, the 
radiation pressure is P r ~ TTF r /c, where L « AnR 2 F r 
in terms of the source radius R and radiation flux F r . 
Thus, since all the balance equations are to first or- 
der linear in F r , it is anticipated that the pressure 
and energy equilibria of the system point to a unique 
value of tt- We explore this possibility in the next 
section. 

3. Pressure Equilibrium For Externally Fed 
Sources 

Let us first suppose that the X-ray source is a 
sphere with Thomson optical depth tt, and that the 
energy is supplied radially by magnetohydrodynamic 
waves. The waves carry an energy density e and prop- 
agate with velocity v a . For definitiveness, we assume 
that these are Alfven waves, in which case the momen- 
tum flux that enters the X-ray source is (1/2) e. The 
magnetic energy of the Alfven waves is in cquiparti- 
tion with the oscillating part of the particle energy 
density, and so we can estimate the gas pressure as 
being of the same order as the ram pressure of the 
oscillating part of the magnetic field, i.e. (1/2) e. Fi- 
nally, we assume that all of the wave energy and mo- 
mentum are absorbed by the source. 

The energy equilibrium for the AR is then given 

by 

F r = ev a , (2) 
whereas in pressure equilibrium 

P r ~ T T F r /c ~ e . (3) 

Dividing the latter equation by the former, one ob- 
tains for the equilibrium Thomson optical depth: 

7"T — ~~~ (4) 

V a 

This value does not depend on luminosity, but it does 
of course depend on the geometry and v a . Nayakshin 
& Melia (1997a) show that the Alfven velocity is close 
to c for typical conditions in an accretion disk if its 
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luminosity is a sizable fraction of the Eddington lu- 
minosity and the magnetic field strength in the AR is 
comparable to the disk thermal energy density. 

Suppose now that the geometry is not perfectly 
spherically symmetric, and that instead the Alfven 
waves can enter the X-ray source through an area A a , 
but the radiation leaves through an area A r ;> A a , 
which is plausibly just the total area of the AR. This 
situation may occur if part of the X-ray source is 
confined by other than the Alfven wave ram pres- 
sure, e.g., by the underlying (non-dynamic) large- 
scale magnetic field (see below). In this case, since 
the energy balance is now F r A r = ev a A ai the equi- 
librium t t is changed to 

To understand the scale represented by the brack- 
eted quantities in this equation, let us consider the 
physical conditions that are likely to be attained dur- 
ing a short-lived and very energetic magnetic flare 
above the standard a-disk. The magnetic field en- 
ergy density is a fraction of the underlying disk en- 
ergy density and the typical size AR a of the flare is 
expected to be of the order of the disk scale height 
(Galeev et al. 1979; Haardt et al. 1994; Nayakshin 
& Melia 1997a). Now, the confinement of the plasma 
inside the flare, and the observed condition 1 > 1, 
require that B 2 /8ir ^> P r 3> P g . Since the magnetic 
stress is much larger than any other stress, the mag- 
netic flux tube adjusts to be in a stress-free vacuum 
configuration. The tube is thick (meaning that its 
cross sectional radius is of the order of its length), 
since the pressure in the disk's atmosphere is insuf- 
ficient to balance the tube's magnetic field pressure. 
This is due to the fact that the magnetic field is pre- 
sumably anchored in the disk's midplane, where the 
pressure is much greater than the atmospheric pres- 
sure. The magnetic waves propagate upwards along 
the magnetic flux tube, while radiation pressure from 
the AR is pushing the gas along the magnetic lines, 
i.e., downwards to the disk. This downward direc- 
tion of the radiation pressure arises naturally in a 
two-phase model (unlike the situation within the ac- 
cretion disk) since here most of the energy is released 
above the disk's atmosphere (see also Nayakshin & 
Melia 1997b) 

With this in mind, we may now describe heuristi- 
cally how the magnetic flare develops and how pres- 
sure equilibrium is established. As is well known 



(Parker 1979; Galeev et al. 1979), magnetic flux tubes 
are buoyant in a stratified atmosphere, and so they 
rise to the surface of the accretion disk. As the tube 
is rising, the particles slide along the magnetic field 
lines downward to the disk in response to gravity. The 
magnetic flux tube becomes more and more particle- 
free, v a is increasing, and so the conditions become 
more and more favorable for the dissipation of mag- 
netic field energy. If magnetohydrodynamic waves are 
generated and propagate up to the top of the flux 
tube, and if the conditions there for effective recon- 
nection are satisfied, the different (incoherent) waves 
reconnect and thus produce highly energetic parti- 
cles. The particles in turn produce X-radiation by 
up-scattering the UV radiation from the disk. Since 
the radiation pressure P r is very much smaller than 
the stress in the underlying magnetic flux tube, we 
may neglect the sideways expansion of the flux tube. 
We need to consider the pressure equilibrium along 
the magnetic field lines, however, since the plasma 
can in principle move freely in that direction. The 
balance of radiation pressure with the magnetic ram 
pressure then sustains the AR optical depth as dis- 
cussed above. Since the flux tube is geometrically 
thick, the corresponding ratio A r /A a is probably of 
order <~ one to a few, and with v a ~ c, we therefore 
expect 

T T ~ 1 - 2 . (6) 

The lowest values of the equilibrium tt can be reached 
due to the fact that A r in this equation is not neces- 
sarily the total area of the source, because some of the 
X-ray flux can be reflected by the underlying disk and 
re-enter the AR. Some of this re-entering flux can be 
parallel to the incoming Alfven energy flux, and thus 
the effective A r is smaller than the full geometrical 
area of the source. Furthermore, we have assumed a 
one-temperature gas, and have neglected the gas pres- 
sure in our calculation. It is possible that the protons 
are much hotter than the electrons, and that they ac- 
count for a sizable fraction of the total pressure in the 
AR, which then leads to a reduction in the value of 
the equilibrium tt as compared to Equation (5). 

In the analysis developed here, there is nothing 
specific to Alfven waves. We could have instead in- 
voked the influx of energy by other waves or even en- 
ergetic particles accelerated by a magnetic reconnec- 
tion process. All that matters is that they have a well- 
defined relationship between their momentum and en- 
ergy densities, and that their propagation speed is 
close to c. Therefore, there may occur a different 
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physical setting where the mechanism for the pres- 
sure equilibrium described here can be at work. It 
seems unavoidable that magnetic fields are involved, 
and that they therefore provide some of the confine- 
ment, e.g., in a direction perpendicular to the mag- 
netic field lines. We encourage future work on the 
range of possible geometries for the ARs. 

4. The Range in tt Permitted by Current 
Observations 

Zdziarski et al. (1996) produced a fit of the average 
Ginga/OSSE spectrum of Seyfert 1 galaxies assuming 
that the active regions form hemispheres above the 
disk. They found that the radial optical depth of the 
hemispheres is tt ~ 1. Here, we will examine whether 
the Seyfert spectrum can be due to a combination of 
spectral components from flares with different tt, but 
the same y-parameter (set arbitrarily at 1.3). The lat- 
ter assumption is introduced to ensure that the X-ray 
spectral index does not vary considerably from flare 
to flare. A constant y-parameter is a natural conse- 
quence of the fixed geometry of the flare, in the sense 
that the cooling of the plasma is fixed by how much 
of the X-ray flux re-enters the emitting region after 
it is reflected from the disk (see Haardt & Maraschi 
1991). 

As an illustration of the method, we first com- 
pute the spectrum from flares with a range of Thom- 
son optical depths assuming that they all have the 
same luminosity. We then convolve these spectra with 
a Gaussian probability distribution that a flare oc- 
curs with tt- The composite spectrum F(E) (in en- 
ergy/sec/keV) is 



F{E) 



cItt exp 



(TT - To) 



F(E,r T ), (7) 



where F(E,tt) is the spectrum from a single flare 
with tt- We take To = 1.14 and adopt several val- 
ues of T a = const to represent the possible spread in 
tt between different flares. The individual spectra 
are computed assuming a slab geometry using an Ed- 
dington frequency-dependent approximation for the 
radiative transfer, using both the isotropic and first 
moments of the exact Klein- Nishina cross section (Na- 
girner & Poutanen 1994). Although this geometry is 
clearly different from that of a realistic flare, our point 
here is to test the possibility of co-adding spectra with 
different tt, in order to see what range in tt may be 
permitted by current observations. We expect that a 



more accurate calculation with the correct geometry 
will yield qualitatively similar limits on tt, though 
the exact values should be inferred using a \ 2 fit to 
the data. 

Figure 1 shows the results of our calculation for 
t = 1.14 and two values of T a : 0.7 and 1.5. The spec- 
trum from a single flare (solid curve) with tt = 1.14 
is also shown for comparison. It can be seen that the 
plot for T a = 0.7 is hardly distinguishable from that 
for tt = 1.14. Moreover, these curves differ the most 
above 100 keV, where the OSSE data typically have 
error bars larger than this difference (see, e.g., Fig. 1 
in Zdziarski et al. 1996). On the basis of this simple 
test, we would expect that Seyfert spectra may be 
comprised of contributions from many ARs encom- 
passing a range (0.5 — 2) of tt- We note, however, 
that a broader range in tt is unlikely because of the 
considerable flattening to the spectrum for T a > 1.5. 

The conclusion that tt is allowed to vary within 
the range of 0.5 — 2 is very important for the magnetic 
flare model, since it is otherwise difficult to see how 
different flares could produce exactly the same tt- It 
may also happen that a flare evolves through many 
phases and that its Thomson optical depth therefore 
varies with time. However, these calculations demon- 
strate that as long as that variation is restricted to the 
range ~ 0.5 — 2, the resulting spectrum is consistent 
with the observations. 

This is not to say that flares with tt greatly in 
excess of 1 are not permitted. According to Equation 
(5), a large optical depth probably corresponds to a 
flare with v a <C c. The magnetic field is therefore 
probably smaller and the energy is transported with a 
smaller v a . That is, magnetic flares with large tt may 
be significantly dimmer than the more energetic ones 
with tt ~ 1, and so they do not contribute noticeably 
to the observed Seyfert spectrum. 

5. Conclusions 

We have considered the consequences of imposing 
a pressure equilibrium on the active regions of Seyfert 
Galaxies, in addition to the more often studied energy 
equilibrium, under the assumption that the emission 
arises within energetic magnetic flares above the sur- 
face of a cold disk. We showed that if the energy 
is supplied to the X-radiating plasma by the influx of 
some energy source with a group velocity <~ c, then tt 
probably falls within the range 1 — 2. The current X- 
ray/7-ray observations are consistent with this range 
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of Thomson optical depths. We conclude that mag- 
netic flares on the surface of the cold disk remain a 
viable explanation for the spectra observed in Seyfert 
Galaxies. Alternative explanations, based on a grav- 
itational confinement of the ARs, appear to be inca- 
pable of accounting for the observed 'universality' in 
the value of tt- 
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Figure 1. The spectrum resulting from co-adding 
the components due to different magnetic flares with 
a Gaussian distribution in the Thomson optical depth 
(see text), centered on ^ = 1.14 with a width T a — 
0.7 (long-dashed curve) and 1.5 (short-dashed curve). 
The spectrum from a single flare with optical depth 
tt = T"o is shown as a solid curve. The curves are 
normalized to the same integrated flux. These spec- 
tra demonstrate that current observations may not be 
able to differentiate between a single-flare spectrum 
and one comprised of many different flares if their 
optical depth is in the range ~ 0.5 — 2. 
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